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In this paper we investigate the non-perturbative generation of an anomalous magnetic moment 
for massless fermions in the presence of an external magnetic field. In the context of massless QED in 
a magnetic field, we prove that the phenomenon of magnetic catalysis of chiral symmetry breaking, 
which has been associated in the literature with dynamical mass generation, is also responsible for 
the generation of a dynamical anomalous magnetic moment. As a consequence, the degenerate 
energy of electrons in Landau levels higher than zero exhibits Zeeman splitting. We explicitly 
report the splitting for the first Landau level and find the non-perturbative Lande ^-factor and 
Bohr magneton. We anticipate that a dynamically generated anomalous magnetic moment will be a 
universal feature of theories with magnetic catalysis. Our findings can be important for condensed 
planar systems as graphene, as well as for highly magnetized dense systems as those forming the 
core of compact stars. 
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I. INTRODUCTION 

The quantum mechanical description of charged fermions in a constant magnetic field has attracted much attention 
since the early developments of Quantum Electrodynamics (QED) The partial Lorentz symmetry breaking 

produced by the magnetic field is reflected in the quantum mechanical properties of the charged fermions which 
behave as free particles along the field, but have quantized momenta Q (characterized by a Landau level /) in the 
transverse direction. At the tree level, the Landau levels I > 1 are degenerate with respect to the spin projections 
along and opposite to the magnetic field. However, this degeneracy is broken by radiative corrections giving rise 
to a spin-field interaction term in the effective action and thus to an anomalous magnetic moment. 

The theory of the electron magnetic moment has historically played an important role in the development of QED. 
As it is known, the electron intrinsic magnetic moment fj, is related to the spin vector s by /x = gfJ-ss, where 
fj-B = eh/2mc is the Bohr magneton, and g is the Lande ^-factor. One of the greatest triumphs of the Dirac theory 
was the prediction of the value g — 2 for the Lande g-factor of the electron in the non-relativistic limit. However, 
this prediction was later challenged by more refined experimental measurements showing a value slightly larger than 
2. The solution of this puzzle came only after Schwinger calculated the first-order radiative correction to fi due to 
the electron-photon interactions [3j. Schwinger's results led to an anomalous magnetic moment with a correction to 
the g- factor of order = a being the fine-structure constant. Higher-order radiative corrections to g have 
subsequently led to a series in powers of a/ir Q that is in excellent agreement with the experiment. 

For massless electrons, however, no anomalous magnetic moment can be found through Schwinger's perturbative 
approach. The problem is that an anomalous magnetic moment would break the chiral symmetry of massless QED, 
but this symmetry is protected against any perturbatively generated breaking term. However, no protection exists 
against non-perturbative breaking of chiral symmetry. A non-perturbative chiral symmetry breaking mechanism is 
known to exist in theories of massless fermions in a constant and uniform magnetic field Q . It is called the magnetic 
catalysis of chiral symmetry breaking (MC\SB). The MC\SB gives rise to a chiral condensate that in turns can 
produce a dynamical fermion mass. Until recently, all the studies of MC\SB focused on the generation of a dynamical 
mass but ignored the possibility of a dynamically generated anomalous magnetic moment. In a recent letter 

we reconsidered the MC\SB in massless QED and showed that in addition to the dynamical mass, the magnetic 
catalysis simultaneously produces a dynamically generated anomalous magnetic moment. Our results gave rise to a 
non-perturbative Bohr magneton proportional to the inverse of the dynamical mass, in analogy to the way in which 
the bare Bohr magneton depends on the bare mass. In the present paper we develop in detail the calculations that 
led to the magnetically catalyzed anomalous magnetic moment found in Ref. [HI]. We also discuss some elucidating 
points about the infrared dynamics of fermions in the lowest Landau level (LLL) and their role in the generation of 
the anomalous magnetic moment for fermions in higher Landau levels. 

The plan of the paper is as follows. In Sec. 2, we present a brief historical account of the appearance of magnetic 
moment contributions in relativistic and non-relativistic theories of charged fermions, including a summary of our 
own findings for the case of massless fermions, which will be developed in detail in the rest of the paper. In Sec. 3, 
we calculate the fermion full propagator in massless QED in a magnetic field, considering the Dirac's structures for 
mass, magnetic moment and wave-function renormalization term. In Sec. 4, an infinite system of coupled Schwinger- 
Dyson (SD) equations for the fermion self-energy is derived within the ladder approximation. This system is then 
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consistently solved in Sec. 5, where explicit solutions of the dynamical quantities are found in the first and in the 
LLL. The dispersion relations for fermions at different LLs in the chiral-condensate phase are calculated in Sec. 6. In 
particular, the Zeeman effect for fermions in the first LL is analyzed and the corresponding non-perturbative Lande's g- 
factor and Bohr magneton are identified. In Sec. 7, we calculate the chiral condensate and find the connection between 
this order parameter and the rest energy of the electrons in the LLL, which is given by the sum of the dynamical 
mass and the magnetic energy due to the anomalous magnetic moment term. This result reflects the infrared origin of 
the condensation phenomenon and underlines the fact that both parameters, the mass and the anomalous magnetic 
moment, can be induced once the chiral symmetry is broken by the condensate. Possible applications of the outcomes 
of this paper for planar condensed matter systems as graphene and for dense astrophysical systems as compact stars 
or magnetars are outlined in Sec. 8. Finally, using a chiral-spin representation for fermions in the LLL, it is shown in 
the Appendix that the dynamics of such fermions reduces to that of free particles in a (l+l)-dimcnsional space with 
an induced rest-energy that depends non-perturbatively on the coupling constant and applied magnetic field. 



II. MAGNETIC MOMENT HISTORICAL REVIEW 



For the sake of understanding we will briefly review in this section the role that spin-field interactions have played 
in non-relativistic and relativistic theories of spin-^ charged particles in a magnetic field. We discuss massive and 
massless QED and explain why in the massless theory we cannot follow the perturbative approach that gave rise to 
an anomalous magnetic moment in the massive situation. It will become clear throughout the paper that a non- 
perturbative method is required in the massless case because the anomalous magnetic moment there will be closely 
connected to the breaking of the original chiral symmetry, hence it has to be dynamically generated through some 
non-perturbative mechanism. 



A. Non-Relativistic Case 



Let us start with the case of a non-relativistic spinless charged particle in the presence of a constant and uniform 
magnetic field H. Assuming a magnetic field pointing along the 23-direction, the energy eigenvalues Es of the 
Schrodinger Hamiltonian 

H s = ^(p-eA) 2 (1) 

are given by @, [l2j 

E s = ^ + ^(n+ 1 -) n = 0,1,2,..., (2) 
2m m 2 

Notice that the particle can freely move along the field direction, while in the plane perpendicular to the field its 
motion is confined to quantum orbits labeled by the discrete numbers n. Without loss of generality, we can assume 
from now on that \eH\ = eH. 

For charged fermions, the Schrodinger Hamiltonian needs to include a spin-field interaction term originally intro- 
duced by Pauli as 

Hp = -^(p-eA) 2 -/x-H, (3) 
2 m 

where 

H = g[i B s (4) 

is the particle intrinsic magnetic moment. Here \xb is the Bohr magneton, s the spin operator given in terms of the 
Pauli matrices as s = er/2, and g the Lande g-factor. For bare electrons g — 2. 
The corresponding energy eigenvalues are 

E P = p- + —(n+h- n B aH 7i = 0,l,2,..., a = ±1 (5) 
2m m 2 

Here a denotes spin projections along (+) and opposite (— ) to the field direction. The last equation can be rewritten 
as 

E P = & + 2Mn+l-l)H=& + 2v B lH (6) 
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where 

l = n+ 1 --^ Z = 0, 1,2, ... (7) 

The number I is known as the Landau level number. Notice the double degeneracy of all the I ^ due to the two spin 
projections contributing to each of them. The LLL (I — 0) is not degenerate because there is only one combination 
of the non-negative integer n and spin projection a that can produce I — 0. As we will see below, the degeneracy of 
/ 7^ can be broken by radiative corrections which give rise to an anomalous magnetic moment term in the electron 
self-energy. 



B. Relativistic Case 



In the relativistic case we should start from the Dirac equation in the presence of an external magnetic field 

(IT^ - m)ijj = (8) 

where the field is introduced through the covariant derivative 

n p - id^ - eA^ (9) 

Assuming again a uniform and constant magnetic field along the X3 direction, and using the gauge A 2 = Hx\, 
Aq — A\ — A3 — 0, Eq. J8|) can be solved for ip. It follows from ([8]) that the energy eigenvalues for the relativistic 
spin-i particle are found [J] to be 



E R = ± s jm?+pl + 2eH{n+^-^) n = J l,2,..., o = ±1 (10) 
with (+) for particles and (— ) for antiparticles. In terms of the Landau levels (JT]) the relativistic energy becomes 



E R = ±^Jm 2 +pl + 2eHl 1 = 0,1,2,... (11) 

Clearly the double spin degeneracy of the I states is also present for relativistic fermions. 

The non-relativistic energy can be easily recuperated by taking the limit 2eH/m 2 <C 1 and •p\jm 2 <S£^ 1 in (fTTjl and 
subtracting the rest energy m, 



E NR = lim E R -m= lim m\/l + -^-H -I - m ~ -P~ + 2a B lH. (12) 

2eH/m 2 ,pl/m 2 ^0 2eff/m 2 ,p§/m 2 - + V m z 171 2m 

Since the spin is automatically incorporated in the relativistic treatment, the non-relativistic limit leads directly 
to the Pauli theory ©. Even more, it naturally gives that a unit of spin angular momentum (tr/2) interacts with 
the magnetic field with a coupling of 2/is, that is, the Dirac theory automatically produces the correct g-factor of 2, 
something that was a puzzle at those times. 



C. Anomalous Magnetic Moment 

Despite Dirac theory's success in predicting g — 2, this result was later challenged by experimental findings of g > 2 
for electrons/positrons. The solution to this puzzle was provided by Schwinger. In a classical paper on the topic 
Q , Schwinger calculated the one- loop contribution to the fermion self-energy in a weak magnetic field that led to an 
anomalous magnetic moment 

H A = (g 1 - 2)fx B s (13) 

Accordingly, the Lande g-factor was modified as 

<?' = 2(1 + ^) (14) 



in good agreement with the experiment. 
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Schwinger's result contained the first order correction in a/%. Subsequent higher-order corrections to g give rise to 
a series in powers of a/n ||. The corrections up to the eighth order has shown an agreement with the experimental 
value that is good to one part in 10 12 [HI]. 

Taking into account Schwinger's anomalous magnetic moment, the Dirac equation acquires an extra structure 
cr^F^ (with = §[7^7*]), so 

(n^y - m + K/J, B HE 3 )ip = (15) 

where K—a/2-K and S3 = 17172 is the spin operator. The corresponding relativistic-particle energy [l4| is 

E l,a = [(m 2 + 2eHl) 1/2 - fi B nH<7} 2 +pl, 1 = 0,1,2,..., a = ±1 (16) 

Hence, once the anomalous magnetic moment contribution is considered the spin degeneracy is removed since the 
energy (|16p explicitly depends on the spin projection a. It is timely to clarify here a mistake appearing in Ref. [l5l |. 
where it was argued that a magnetic moment term (oy^-F^) could not be present in the self-energy because it would 
give rise to an energy that would depend on the orientation of the magnetic field. It is obvious from (fl6|) . that no 
matter what the direction of the magnetic field is, the energy of the particle with spin oriented along or opposite 
to the field will not change. Notice that a 180° rotation equally affects the magnetic field and the spin, leaving the 
magnetic moment contribution to the energy — fj, ■ H unchanged. 

As was done in Eqs. (H2), we can take the non-relativistic limit of (jXSJ) to find 

E NR = ^ + 2(n+ \)n B H - 2(1 + n)^H = 

= + 2(n + hfi B H - g'^H = ^- + 2lfi B H - n^ B aH, (17) 
2m 2 2 2m 

This result shows that the Lande g-factor depends on the fine-structure constant as pointed out in (|14p . and that the 
anomalous magnetic moment breaks the spin degeneracy of all LLs with I > 1, hence producing the following energy 
splitting of the levels 

AE t = 2^ B H = (g' - 2)fi B H (18) 

The case with strong fields, i.e. fields about B c ~ 10 13 G, was considered in Refs.fl6l|-[l7l]. In those works, the authors 
treated the field exactly, (i.e. without expanding in powers of H), in the fermion one- loop self energy but kept a 
perturbative treatment in the coupling constant. In that approximation the energy splitting due to the anomalous 
magnetic moment no longer changes linearly with H, and besides it depends on I [l7| . 



D. Induced Magnetic Moment for Massless Fermions 

As previously shown, the electron magnetic moment is related through the Bohr magneton [i B to its charge-to- 
mass ratio. Hence, the origin of the extra contribution due to the anomalous magnetic moment appearing in quantum 
field theory can be understood from the fact that there the electron is continuously self-interacting through radiative 
interactions. Thus, part of the electron energy, and consequently of the electron mass, will be transferred to the 
created photon cloud. Therefore, as a consequence of the decrease of the electron's mass, the corresponding magnetic 
moment will be strengthened. 

On the other hand, as we have already stressed in the Introduction, if one is interested in exploring the appearance 
of an anomalous magnetic moment in a theory of massless fermions, a perturbative, "a-la-Schwinger" approach is not 
possible anymore. Since an anomalous magnetic moment would break the chiral symmetry of the massless theory, it 
can only be generated via non-perturbative effects. 

Henceforth, we are going to explore the dynamical generation of an anomalous magnetic moment in the context of 
massless 4-dimensional QED in the presence of a uniform and constant magnetic field. As we are going to show in 
the next sections, the formation of a chiral condensate through the MC\SB mechanism is responsible not just for 
the dynamical generation of a mass, but also for the appearance of a dynamical magnetic moment. 

Physically it is easy to understand the origin of the new dynamical quantity. The chiral condensate carries non-zero 
magnetic moment, since the particles forming the condensate have opposite spins and opposite charges. Therefore, 
chiral condensation will inexorably provide the quasiparticles with both a dynamical mass and a dynamical magnetic 
moment. Symmetry arguments can help us also to better understand this phenomenon. A magnetic moment term 
does not break any additional symmetry that has not already been broken by a mass term. Hence, once MC\SB 
occurs, there is no reason why only one of these parameters should be different from zero. 
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We will show below that a very important consequence of the dynamically generated magnetic moment is a splitting 
in the electron energy spectrum that can be interpreted as a non-perturbative Zeeman effect. In the LLL, since only 
electrons with one spin projection are allowed, there is no energy degeneracy and therefore no splitting can occur. 
However, for electrons in higher LLs the energy of the degenerated spin states is splitted by the interaction of the 
induced magnetic moment with the applied field. The corresponding energy splitting can be conveniently written in 
the well known form of the Zeeman energy splitting for the two spin projections as 

AE = 2np, B H (19) 

where k and pb are the non-perturbative Lande <?-factor and Bohr magneton respectively. For electrons in the first 
LL they are given by [ll[ 

* = e" 2 ^, MB = ^ (20) 

Worth to highlight in the above results are the non-perturbative dependence on the coupling constant a in the Lande 
g- factor and the Bohr magneton's dependence on the dynamically induced electron mass M 1 . 

III. ELECTRON FULL PROPAGATOR IN MOMENTUM SPACE 

In this Section we are going to obtain the fermion's full propagator in QED with masslcss bare fermions in the 
presence of a constant and uniform magnetic field. The full propagator obeys the following equation 

p^" - £(>, x')]G{x, x') = 8 A {x - x') (21) 
The structure of the the self-energy T,(x,x') [l8[ 

£(z, x') = (Z||nJI 7 f + Z x Uj^ + M + ^F^a^)5\x - x') (22) 

contains the wave function's renormalization coefficients Zn and Z_\_, as well as mass M and anomalous magnetic 
moment T terms, all of which have to be determined self-consistently as the solutions of the SD equations of the 
theory. 

In (f2"2"| F^ v — F^ /H denotes the normalized electromagnetic strength tensor, with H the field strength. The 
external magnetic field breaks the rotational symmetry of the theory, hence separating between longitudinal p\\ ■ 7" — 
V V K P F* W 1^ (for (i,v = 0, 3), and transverse p ± • 7 X = (for //, v = 1, 2), modes. F* v = ^ n e^ pX FP x is 

the dual of the normalized electromagnetic strength tensor F} t „. 

The transformation to momentum space of (f2"2")) can be done by using the so-called Ritus' method, originally 
developed for fermions in [lj| and later extended to vector fields in [2(|. In Ritus' approach, the transformation to 
momentum space is carried out using the eigenfunctions E l p (x) of the asymptotic states of the charged fermions in a 
uniform magnetic field 

E l p (x) = E+(x)A(+) + E-(x)A(-) (23) 

where 

A(±) = (24) 

are the spin up (+) and down (— ) projectors, and Ep^ (x) are the corresponding eigenfunctions 

E+(x) = N l e-^ PoX ° +P2x2+P3X ^ Di(p), 

E-{x) = Ni^e-^^+P^+P^Dt^ip) (25) 

with normalization constant Ni = (AneH) 1 ^ J 'y/l\, and Di(p) denoting the parabolic cylinder functions of argument 
p = \/2eH(xi — p2ls.II), and index given by the Landau level numbers I = 0, 1, 2, .... 
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The E l functions satisfy the orthogonality condition [l5| 

J d 4 xE l p (x)E 1 ;, (x) = (27r) 4 ?W (p - p')U(l) , (26) 

with^=7 (4)t 7 o, 

6^ (p - p') = S u 'S(p - P ' )S(p2 - p' 2 )5(p 3 - p' 3 ), (27) 

and 

11(1) = A(+)S 10 + 1(1 - S 10 ). (28) 
The spin structure of the E p functions (|2"3")) is essential to satisfy the eigenvalue equations 

{n- 1 )E l p {x) = E l v {x){ 1 -p) , (29) 

and 

(^njl 7 f + Z x n^)E l p (x) = E l p (x)(Z {{ p^l + Z x p^) , (30) 

with^ = (p°,0,-V2eHl,p 3 ), thusp^ = (0, 0, -\f2eHl, 0) and = (p°,0,0,p 3 ). 

The relations (l29 |) - ([30|) and the orthogonality condition (|26| . facilitate the diagonalization of the fermion self energy 
£(x,a/) in momentum space 

E(p,j/) = J d i xd i yE l p (xMx,y)E p :(y) = (2^6^ (p - p')Il(l)X l (p) (31) 

with 

E< (p) = Zf pjf 7 Jl + Zi^ 7 ^ + M'l + ^7 V (32) 
Using the spin projectors (|24p and introducing the longitudinal and transverse projectors 

Aj= = i(l±2^p!), A± = i(l±z 7 2 ). (33) 

the function E'(p) can be rewritten in the following form, 

E<(p) - Zjj(A+ - Ap|p„| + iZi(Aj - A+)|P_l| + (M ; + T ( )A(+) + (M« - T')A(-) (34) 

It is clear from (|3ip that when the transformation is correctly used (i.e. taking into account the projector H(l) 
in the orthogonal condition), the separation between the LLL and the rest of the levels is automatically produced. 
Considering I = in (|34|) . and using that p±(l = 0) = and A(+)A(— ) = 0, the spinorial structure in the RHS of 
Eq. (f3"Tj) reduces to 

n(0)E°(p) = [Z°(A+ - Ap|Pu| + (M° + T°)A(+)]A(+) (35) 

While at I ^ 0, it is given by 

11(1 + 0)E'(p) = f}{p) = Z|(Ajj--Ajj-)|p|||+iZi(Aj- A+)|JJJ + (M l +T l )A(+) + (M l - T l )A(-) (36) 

From a physical point of view the results (|35j) and (|36p are simply reflecting the fact that there is only one spin 
projection in the LLL , and hence the self-energy E (p) only contains the spin projector A(+). Since the remaining 
LL's contain two spin projections, E(p) depends on the two projectors A(+) and A(— ). 

With the help of Eqs. P6"|) . (|!?9")) - (f3"l"j) . it is straightforward to show that the inverse of the full fermion propagator 
in momentum space is given by 

GrW) = / rfV^Wp.T-^.^fe) = {2^Xp-p'm)[p-l-^{p)\ (37) 
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The full propagator G l (p,p') must satisfy 

^0^G^(p,pnG r \p»,pi) = (2n)^(p-p')Il(l) (38) 

where XpDj^ = Ez dpo ^ P3 . It is easy to see that §8$) is indeed satisfied by 

G l (p, p') = (2tt) 4 ?( 4 ) (p - p')Il(l)G l m (39) 

with G l (p) formally given by 

G l ® = ~ ^-r (40) 

To find the explicit form of G l (p) we have to solve 

G\p)Gi 1 (p) = Gi\p)G l (p) = I (41) 

where 

Gi\p) = P ■ 7 - = 7 • V 1 - M l I - THw (42) 



and V M = ((1 - Z|)po,0, (1 - Z^VMTl, (1 - Z\ { )pa). 
One can show that the matrix function 



with 



and 



G l m = - % u ( 43 ) 

detG, 



A = 7 1 Gr- 1 (p) 7 i, B = 75Gr(p)^7s (44) 



det G ; (p) = \jdct[Gi l {p)AB] 

= -{[M l + (Vjf - T') + V[][M ; - (Vjf - T') + V_[] 

+ [Af ( + (Vjf - T*) - Vj][M' - (Vjf - T l ) - V\_]} 

x {[M l + (Vjf + T 1 ) + Vjj[M* - (Vjf + T z ) + V|] 

+ [M ( + (Vjf + T l ) - V[][M l - (Vjf + T ; ) - Vj]}, (45) 

satisfies the condition (F4"T]) . 

Using (|4"3")) and working in the basis of projectors (j2~4")l and (|3"3")l . the matrix function G l defining the full fermion 
propagator (jST))) can be written as 



with notation 



A^ Z (T,^|) = T l - M l - Vjf 

L>'(T) = (Af') 2 - (Vjf - T 1 ) 2 + (Vj) 2 

Vjf = 

V[ = (1 - Zi)|p x | = (1 - Z{)VteM. (47) 
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In the LLL, = 0, thus the LLL full propagator becomes 

+ vf-(/f°-r°) A( - )A ii - vf + (^-r°) A( " )A ii (48) 

IV. SCHWINGER-DYSON EQUATION FOR THE FERMION SELF-ENERGY 

To explore the dynamical generation of a magnetic moment in massless QED, we can start from the SD equation for 
the fermion self-energy in the presence of a constant magnetic field. We will work in the quenched-ladder approximation 
where 

E(x, x') = -ieVG(i, x')j"D^(x - x'). (49) 

Here, E(x, x') is the fermion self-energy operator (f2"2")l . D^ v {x — x') is the bare photon propagator, and G(x, x') is the 
full fermion propagator depending on the dynamically induced quantities and the magnetic field. 
Equation (|49[) can be transformed to momentum space by using the E p functions as 



d 4 xd 4 x'E l p (x)E(x,x')E l p ,(x') = -ie 2 j d 4 xd 4 x' ' E p {x)^ 
(^^^.(x)n(r0G r (r)4'( a;, ))7 !/ <'(^)^^-ar') 



(50) 



where 



, d A n e -iq-( x - x ') a a 

D K»( X -*')= <o V 2 ■ ^ - 0- ~ 0^)> (51) 

' (27r) 4 q A — ie g z 



with £ the gauge fixing parameter, and we used that 



G (^') = ^^W n 0")G'>')4 : M (52) 



At this point it is convenient to consider the integrals 
-/ 



d 4 xE p {x) 1 tl E p l{x)e- lq - x = (2-k) 4 S (3) (p^ + q ^ p)e~ l(!l{p2 " +P2)/2eH e~^ /2 

x 7^=^4.(^(^(0"), (53) 
— ' vn!n ! 



and 



dV^^OT'^'/CxOe^- 35 ' = (2tt) 4 ^(p" +q- 39 ')e l9l(p " 2+pi)/2eff e-^ /2 



x 



£ ^i==e l (""-"')^J^ n ,(^)A(a") 7 ^A( ( T'), (54) 



with n = n(Z, <x), n" = n(7" , a"), n' = n(Z ; , a'), and n" = n(Z" , a"), defined according to 



„(/, tr ) = i + |-i / = 0,1,2,..., <7 = ±1. (55) 



The notation in (|53j) and (J54j) included the use of polar coordinates for the transverse g-momentum q± = \Jq 2 + qf , 
95 = arctan($2/<2i); normalized quantities = Q fJ _/\ /, 2eH; the tri-delta function 

5^ (p» + g - p) = <5(p " + 90 - Po)% 2 " + 92 - p 2 )<$(p 3 " + 93 - p 3 ); (56) 
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and 



min(n.n" ) . , . 

"1. T) , I 



J n „»(«i) = £ -fi ryr^ ^^]" + ""- 2m . (57) 

z — ' m!(n — m)!(n — my. 

Doing the integrals in x and x' in (|50|) with the help of (|53[) and (|54[) . integrating in p v , and using the Feynman gauge 
(£ = 1), one finds 

cfiq y y gi(n— n"+n"— n')v g-gj 



Vn!n'!n"!n"! q 2 

x4„-.(gi)Jrr l '(9i)A(^A(^)n(r)G l "(F9)A(^")7' J A(^) I (58) 



where p — q = (p° — <?°, 0, —y/2eHl" , p 3 — g 3 ) and [er] means summing over er, a', a" , a" . The appearance of the II(Z) 
factors in both sides of the equation ensures the correct counting of only one spin projection for the fermions at the 
LLL. 

Due to the negative exponential e~ q± the main contribution to (|58|) will come from the smallest values of q±. This 
allows one to keep in ([55)1 only the terms with the smallest power of q± in Jn n n'(q±) (see Ref. [8] for details). Hence 

[max(n',n")]\ i „i 

Jn„»(2-L) -> _^ [Zgx] 1 1 -> ™!<W (59) 



and we obtain 



£ , (»n(i)«« , =-fc 9 (2efl)/^x:x; e " sl 



r [<t] 

x^„„^ ?r .„,A( ( T)7 /i A( ( T")n(r)G r (F 7 g)A(a") 7 ^A( ( T'), (60) 

Taking into account that 

Sn,n" = Sl t l"5 attT " + 5j+<r,i" £-<r,<7" (61) 

together with the relations 

A(±)<y£ = <A(T), A(±) 7 j! = 7 J!A(±) (62) 

A(±)A(±) = A(±), A(±)A( T ) = 0, 7 ,YtJ=0, (63) 
we can do the sums in [a] and I" in (foT))) to arrive at the SD equation 

E<(p)n(0 = -ie 2 (2eff)n(0 | ^/^ [ 7 llff(p = g) 7 ll 

+A(+) 7 ^G i+1 (F^)7^A(+) + A(-) 7 ^G i - 1 (F^)7 M ± A(-)] (64) 
If the external fermion is in the LLL (^ = 0), Eq. reduces to 

S°(p)A(+) = - i e 2 (2eH)A(+) J ^^i^^bjl + ^G 1 fr=q)^ (65) 

with S°(p)A(+) given in Eq. ([331) . While if the external fermion is in any higher LL (I ^ 0), Eq. (|G4[) becomes 

£^°>(p) = -ze 2 (2etf) / ^^[^(F^b)! + A(+) 7 ^G i+1 (F^A(+) 

+A(-) 7p i G 1 - 1 (F^)7, 1 A(-)] (66) 

with £^°)(p) given in Eq. flM]). 
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Since the equation for a given Landau level I involves dynamical parameters that depend on I, I — 1 and I + 1, 
the SD equations for all the LL's form a system of infinite coupled equations. Fortunately in the infrared region, 
the leading contribution to each equation will come from the propagators with the lower LL's, since the magnetic 
field term (~ IB) in the denominator of the fermion propagator for I ^ acts as a suppressing factor. Using this 
approximation, one can find a consistent solution at each level. It is also convenient to notice that, since the solution 
for the first LL depends on the one for the LLL; the solution for the second LL depends on the one for the first 
LL, and so on, the solutions for M l and T l can all ultimately be expressed as a function of the LLL solution. This 
indicates that the physical origin of all the dynamical quantities is actually due to the infrared dynamics taking place 
at the LLL. 



V. INDUCED ELECTRON MASS AND ANOMALOUS MAGNETIC MOMENT 

A. Solution of the SD Equation in the LLL 

Let's work in the LLL, that means, p± = 0, and consequently the LHS of is given by ([55)1 . In the infrared region, 
the leading contribution to the RHS of Eq. (|65[) comes from the term with no magnetic field in the denominator, that 
is, the term with the LLL propagator G° (p — q). Hence, in the leading approximation Eq. (|65[) is given by 

(Af + r°)A(+) + Z°A(+)(A+-A-)||5 || |^-ze 2 (2ei/)A(+) J ^ e g x 7 jlffl(p = g) 7 ll (67) 



II 

Using (|6"2 ]) -(fgg ]) together with 



l AH = 1. (68) 



assuming <C 1, and doing a Wick's rotation to Euclidean space, we obtain 

(M° + T°)A(+) + Z°A(+)(A+ - Ap|p|| | = e 2 (2eH) 
K . , [ d 4 q (M° + T°) 

From Eq.fjnU) it is clear that Z9 = 0, which corroborates the assumption we did before. Taking the infrared limit 
(pu ~ 0), and assuming that M° + T° is independent of the parallel momentum, we arrive at 

If M° + T° is replaced by md yn , Eq. ([70)) becomes identical to the SD equation obtained in the phenomenon of 
magnetic catalysis of chiral symmetry breaking. Thus, the solution of (|70[) is given by 



M° + T° ~ V2e#e~V « (71) 

As in this solution is obtained considering that M° + T° does not depend on the parallel momentum, an 

assumption consistent within the ladder approximation [2~H . As proved in [22| , when the polarization effect is included 
in the gap equation through the improved-ladder approximation, the solution for rrid yn is of the same form as (|71[) . 
but with the replacement y/n/a — > 7r/alog(7r/o;) in the exponent. Since the inclusion of the magnetic moment in the 
LLL SD equation merely implies the replacement md yn —> M° + T°, it is clear that a similar effect will occur in the 
solution (fTTj) . However, this effect will not qualitatively change the nature of our findings. 

As the fermions in the LLL have only one spin orientation it is not possible to find M° and T° independently (see 
that the combination M° — T° is absent from the LHS of the SD equation (|6T|) . as well as from the RHS, because the 
spin projector A(+) ensures that only the terms containing M° + T° in |48|) contribute to the RHS of (|67|)). Thus, 
the SD equation determines the induced LLL rest-energy 

E a = M° + T°. (72) 



which has contributions from the dynamical mass and from the magnetic energy related to the interaction between 
the magnetic field and the dynamically induced magnetic moment. 
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B. Solution of the SD Equation in the First-LL 

The SD equation for a fermion in the first-LL is 

&®=-ie 2 (2eH) J ^ e ^ [ 7 jlffl(p^g) 7 ll + 
+A{+)^G 2 (J— q )^A(+) + A(-) 7 ^G°(F^)7^A(-)] (73) 

Here again the leading contribution comes from the term with no magnetic field in the denominator. Hence to find 
the leading contribution we just need to keep the term depending on G°(p — q) in the RHS. Using (|62|) -([63 |) . together 
with 

7^A± 7 ^ = 2AJ, (74) 
and working as before in the infrared limit, the SD equation after Wick's rotation becomes 

Zl 72 (2 ej ff) + (M 1 +T 1 )A(+) + (M 1 -T 1 )A(-) = e 2 (4 e i/)A(-) J (2 ^ 4 ^ (£0)2 + g 2 (75) 

Therefore, 

M 1 + T 1 = 0, Z\ = (76) 

and 

Mi-rW^/l^g* (77, 



Taking into account ([70)) in (|77j) we get 



M i _ T i = E o ( 7& j 



Finally, from (TFBJ and J78J), it results 



M 1 = -T 1 = ^E a = v / eH/2er^, 



(79) 



The solution |79|) corroborates the relevance of the LLL dynamics (both M 1 and T 1 are determined by E°) in the 
generation of the dynamical mass and magnetic moment of the fermions in the first LL. Given that the magnitude of 
the magnetic moment for the electrons in the first LL is determined by the dynamically generated rest-energy of the 
electrons in the LLL, any modification of the theory producing an increase of E° will, in turn, lead to an increase in 
the magnitude of T l . From the experience with the MC\SB phenomenon, such modifications could be for example, 
lowering the space dimension (23j , introducing scalar- fermion interactions l2l"| , or considering a non-zero bare mass 

HI- 

For the remaining LL's the procedure is similar. For example, the leading term in the second LL will be given by 
the G 1 (p — q) contribution which in turn depends on Eq through the values found for M 1 and T . Therefore, the 
values of M l and T l for higher LL's will depend on E° through the found values of the previous LL's. This fact shows 
that the infra-red dynamics of the electrons in the LLL is the dominant one. 

VI. DISPERSION RELATIONS AND ZEEMAN SPLITTING 

A. Fermions in the LLL 

In the Appendix we showed that the Dirac equation for fermions in the LLL can be written as 

[p. 7 -£°]V> iM = 0, (80) 

where V'lll is a spin-up two-component spinor. Eq. (JSUJ coincides with that of the free (1+1)-D Thirring model [25| . 
with the (1 + 1) — D gamma matrices 70 = oi, 71 = — i(X2, defined in terms of the Pauli matrices ai, and = (po,P3)- 
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The dispersion relation of the LLL fermions obtained from (|80|) is 

Po = ±y/p> + (E°)2, (81) 

Thus, the effect of a dynamical magnetic moment is irrelevant for the LLL fermions, since it just redefines their rest 
energy through the replacement mdyn ~ ¥ M° + T°. This is physically natural, since the fermions in the LLL can 
only have one spin projection, so for them there is no spin degeneracy and hence, no possible energy splitting due to 
the magnetic moment. We shall see below that the dynamical anomalous magnetic moment turns out to be really 
relevant for fermions in higher LL's. 



B. Fermions in the First-LL 



Let us find now the dispersion relations for fermions in higher LL's, taking into account the dynamically induced 
quantities. Starting from the modified field equation in the presence of the magnetic field, 



[p-7-M ; J-iTV7 2 ]^ =0, (82) 

i 
II 



where we neglected the coefficients Z in the terms (1 — Zi) and (1 — Z l ± ), since it is expected that the Z's are much 



smaller than one; the dispersion relations are found from 

det [p ■ 7 - M l I - irV 7 2 ] = [{M l f - feu - T l f+pl][(M l f - (p„ + T l f +p\] = 0. (83) 

yielding- 



ly =^+[V(A^) 2 +2effl±T<] 2 , (84) 

and thus showing that the induced magnetic moment breaks the energy degeneracy between the spin states in the 
same LL (see the double sign in front of T l ) . 

In particular for the first-LL, and taking into account that M 1 /y / 2eH , T 1 j \/2eH <C 1, the leading contribution to 
the energy becomes 



pi ~ p\ + 2eH + (A/ 1 ) 2 + (T 1 ) 2 ± 2T 1 v / 2iff, (85) 

Working in the infrared region (p\/2eH <§; 1) the expression (|85|) can be approximated as 

r , p\ (M 1 ) 2 + (T 1 ) 2 tl , s 

2V2^H 2V2UI 

As a consequence, the energy splitting for the fermions in the first LL is 



AE = I2T 1 ! = 2v / aff/2Vv' r / a (87) 

One can rewrite the last expression in the usual form of the Zeeman energy splitting for the two spin projections 
already introduced in Sec. [ill in Eq. (|19|) as AE = gjlsH with g and JIb representing the non-perturbative Lande 
g-factor and Bohr magneton respectively given by g = 2e _2 V / ^/™, J1 B = T^gr- 



VII. CHIRAL CONDENSATE 



Now we are going to find the relation between the chiral condensate that is, the order parameter of the 

magnetically catalyzed chiral symmetry breaking, and the LLL fermion rest-energy, which as seen before depends on 
the sum of the LLL mass and magnetic moment energy contribution. 

We start from the definition of the chiral condensate 

=iTr[G(x,x)/V], (88) 

with G(x, x) the full fermion propagator given in (|52|) , and V the system volume. As discussed before, the phenomenon 
of magnetic catalysis of chiral symmetry breaking is physically due to the infrared dynamics of the fermions in the 
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LLL. Thus, the leading contribution to the condensate (|88|) comes from the LLL fermions. It is convenient to recall 
here that the space-dependent part of the LLL fermion wave function is given by [26j 

(Xj-X,.) 2 

*(x) ~ S x ° P0+x2p2+xip3) e "b , (89) 

where x c = P2^b is the coordinate of the center of the Landau orbits and Ib = \j\feB is the magnetic length. It is 
clear from (189|) that a particle in the LLL can be localized along the x-i and X3 directions up to infinite, but along 
the x\ direction it is confined within a magnetic length due to the Gaussian function with width Ib (note that for 
(|89[) the standard deviation from the particle position average value (x\) is er = Ib)- This implies that in the LLL the 
volume V in (|88p is given by V = L0L2L3IB, with Li = J_™ dxi. 

Using (|52[) in (|88[) and keeping only the leading LLL contribution in the sum, we obtain after integrating in x's, 

.L L 2 L 3 ^ f f dp dp 2 dp 3 sP?o/-\ 



= i-^-tr{ J ™*y* A(+)G°(jp)} (90) 

with G (p) given in (|48|) . and tr denoting the remaining spinorial trace. In the above result we used the orthogonality 
of the parabolic cylinder functions 

f 00 

dpA(p)A'(p) = v^FlW H / (91) 



The trace operation reduces the previous expression to 

^ l l B J (2tt)3 OT + bnP 

Transforming to Euclidean space, taking polar coordinates for the parallel momenta, and using a momentum cut-off 
defined by the magnetic scale 1/Ib, which is dominant in the infra-red region, we have 



After integrating in the momenta we obtain 



m ^-^ E »H^L) (94) 



This result shows that the role played by the dynamical mass in previous works on magnetic catalysis, on which the 
induced anomalous magnetic moment was ignored Q, is now played by E°. Considering the magnetic- field dependence 
of E° given in (fTTj) , we can rewrite (|94|) as 

E° = IW) (95) 

m2 — Wir/a 

It shows that the induced rest-energy of the electrons in the LLL is proportional to the condensate. This reflects the 
fact that the same order parameter produces the induction of the two quantities contributing to E°: the dynamical 
mass and anomalous magnetic moment. This also confirms that the dynamical mass and magnetic moment have a 
common physical origin, as they are related to one single order parameter, the chiral condensate (^^). 

This validates the fact that since the induction of the term a tlv F ,J ' v in the electron self-energy does not produce any 
additional symmetry breaking to that already created by the mass term, then once the chiral symmetry is broken by 
the condensate (^^'), both dynamical parameters M° and T° are generated. 



VIII. CONCLUDING REMARKS 



In this paper we have presented a fresh view of the phenomenon of magnetic catalysis of chiral symmetry breaking 
in massless QED. Our results show that in this phenomenon the dynamical generation of an additional parameter: 
the fermion anomalous magnetic moment, has to be considered in equal footing to that already taken into account 
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in previous works on this topic: the fermion mass. The rationale for the increase in the number of the induced 
parameters is easy to understand on symmetry arguments. The reason is that once the chiral symmetry is broken by 
the magnetically catalyzed condensation of fermion/anti- fermion pairs, all the physical quantities in the action that 
were forbidden by the chiral symmetry of the bare theory can be now dynamically induced. Therefore, the chiral 
condensate induces dynamical mass and anomalous magnetic moment, since in the massless theory their generation 
was only protected by that same symmetry. Another important outcome of this paper is the lifting of the spin 
degeneracy of all the LLs different from zero. The degeneracy is eliminated due to the Zeeman splitting produced 
by the dynamical anomalous magnetic moment. This energy splitting is given in terms of a non-perturbative Lande 
g-factor and a Bohr magneton that depends on the dynamical mass. 

We call attention to a fact that has not been emphasized in previous works on magnetic catalysis, but that is really 
essential to prove the results here obtained. We refer to the dependence of the dynamical quantities on the Landau 
levels. There is no reason to assume a priori that the dynamical parameters have to be the same for all the LLs. One 
can see that, as soon as one is interested in exploring the SD equations for higher LLs. It becomes clear that there is 
no consistent solution of the infinite system of coupled equations, unless one takes into consideration that the mass 
and the magnetic moment depend on the LLs. This is evident already from Eq. (|75p . 

We anticipate that the dynamical generation of the magnetic moment will be a universal feature of theories with 
magnetic catalysis, in the same way that occurs with the generation of the dynamical mass. Moreover, we expect 
that an anomalous magnetic moment will be also dynamically generated in theories with fcrmion/fcrmion condensates 
in the presence of a magnetic field, as long as the symmetry broken by the condensate coincides with the one that 
would be explicitly broken by a magnetic moment term in the action. These considerations point to some potential 
applications of our findings in two different areas: condensed matter and the ultra dense matter existing in the core 
of compact stars. 

In condensed matter, the most plausible application at present could be in the physics of graphene. It is known 
that the 2-dimensional crystalline form of carbon known as graphene [27j has charge carriers that behave as massless 
Dirac electrons. In particular, a phenomenon where the dynamically induced Zeeman effect we have found here 
could shed some new light is the lifting of the fourfold degeneracy of the I — LL, and twofold degeneracy of the 
/ = 1 LL in the recently found quantum Hall states corresponding to filling factors v = 0, ±1,±4 under strong 
magnetic fields [28] . An attempt to explain the observed lifting of the LLs degeneracy was carried out in [2!| with the 
help of a 2+1-dimensional four-fermion-interaction model of Dirac quasiparticles with chemical potentials interpreted 
as Quantum Hall Ferromagnetism (QFH) order parameters, and dynamical masses related to the phenomenon of 
MCxSB. Nevertheless, we should call attention to two shortcomings of Ref. [2{|. First, the order parameter Jl s was 
introduced in (29| as a chemical potential and interpreted as a QFH order parameter, implying that its physical origin 
was considered unrelated to the phenomenon of magnetic catalysis. However, it is not hard to understand, based on 
the results of the present paper, that /i s should have been identified as a 2+1 dynamical anomalous magnetic moment, 
thus similar to the one we have found in 3+1 massless QED. In other words, ju s should have been connected to the 
MCxSB phenomenon. Second, the LLL results reported in the second paper of [29[ are not self-consistent, since 
they were found taking into account the single (pseudo)spin contribution in the LLL on the RHS of the gap equation 
(A31), but ignoring it in the self-energy operator appearing in the LHS of the same equation. As a consequence, 
separated values for the LLL parameters analogous to our M° and T° were obtained. 

The other possible application in the realm of ultra-dense matter relies on the phenomenon of color supercon- 
ductivity. An important aspect of color superconductivity is its magnetic properties [30l]-[3"3j|. In spin-zero color 
superconductivity, although the color condensate has non-zero electric charge, it is neutral with respect to a modified 
electromagnetism. This so-called rotated electromagnetism is present because a linear combination of the photon 
and the eight gluon remains massless, hence giving rise, in both the 2SC and CFL phases, to a long-range "rotated- 
electromagnetic" field [30} . The long-range field can propagate in the color superconductor implying that there is no 
Meissner effect for the rotated component of an external magnetic field applied to the color superconductor. Even 
though the quark-quark condensate is neutral with respect to the rotated charge, an applied magnetic field can inter- 
act with the quarks of a pair formed by Q-charged quarks of opposite sign and, moreover, for large magnetic fields this 
interaction can reinforce these pairs [32j . In this sense, the " rotated-electromagnetism" in the color-superconductor 
has some resemblance with the chiral condensate in a theory with magnetic catalysis. It is natural to expect then 
that a dynamically magnetic moment can also be induced in a color superconductor under an applied magnetic field. 
Since, on the other hand, the Meissner instabilities that appear in some density regions of the color superconductor 
can be removed by the induction of a magnetic field [331 ] , it will be interesting to investigate what could be the role 
in this process of a dynamically induced magnetic moment. 

Acknowledgments 

This work was supported in part by the Office of Nuclear Physics of the Department of Energy under contract 
DE-FG02-09ER41599. 



15 



APPENDIX A: THE LLL LAGRANGIAN IN THE CHIRAL-CONDENSATE PHASE 

We are interested in obtaining the Dirac equation in momentum space for fermions in the LLL of the chiral- 
condensate phase (|5H)) . With this goal in mind, we should start from the Dirac Lagrangian in the presence of a 
constant and uniform magnetic field that includes the self-energy corrections 



£ = J d 4 ^(x)(n M7 f -£(z))V(z) (Al) 
Using Ritus' transformation to momentum space for the wave functions 

and taking into account (|26l) . (|29|) and (j3Tj) . we obtain 



1=1 

with 



Considering the Dirac matrices in the chiral representation 



where a % are the Pauli matrices, we can express the spin projectors (|24f in terms of a % as 

a 

with 



and introduce the chiral projection operators 



(A2) 



c = ^j^Mpmmrp, - v\p)]Mp) (as) 

The factor 11(7), given in Eq. (128 j) . separates the LLL Lagrangian Co from the rest 

c = CQ+ jr[ d i^i^ {p) [7% _ mi {p) (A4) 

J 2-7T 



Co - I ^^WAt+Jb 1 - PI, - ^(P)]MP) (A5) 



v-^iVo 1 ). M 1 ^), (A6) 



A(±) = r „± ) - (AT) 



^^(lia 3 ), (A8) 



The projectors (|A7|1 and (IA9p satisfy the commutation relations 

[A(±),L] = [A(±),J2]=0 (A10) 
Introducing now the chiral-spin representation for the Dirac spinor 

^ = RA{+)if>, ^^MH^, V>£ +) = LA(+)</>, Vi -) = ^A(-^ (All) 
V4 +) =^LA(+), ^'=#AH, ?i +) -?i?A(+), ^ _) =^A(-) (A12) 
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and using that 



we obtain 



vb = v4 +) + A ] + A +) + i> [ L ] (A13) 



Co = [ *^*»fi&>(p) 7 « . p ^(p) -^g(p)^W(p) + V&W 1 -Pii^W -Sf(p)^(p)] (A14) 
J (27r) 

where we used A(+)S°(p) = £°A(+) since Z9 = (sec Eq. (|55j)). From (|A14j) we see that the LLL only gets 

contribution from the wave functions of the spin up-states. Introducing ipj — (ipi, ip2i ip3, 1P4) for the LLL four spinors 
we can rewrite Eq. (|A14|) as 

^-/^w.w^ <*+»>)(*). (a 15 , 

Notice that in the LLL the fermion spinor reduces to a bispinor which corresponds to the two chiralities of the spin 
up state in this case. 

At this point it is more convenient to work with the (1+1)-D gamma matrices 

t°=oi= (; ;) , ?=-ia 2 = (; - 1 ) (ai 6 ) 

which satisfy the algebra 

= + e^Y, (A17) 

7 M 7 5 = ~^ v lu (A18) 
with 7 5 = 7°7 1 and the (1+1)-D metric and the totally antisymmetric tensor given respectively by 

i A ,u ( o i 



-1 I ' " I -1 



(A19) 



1>lll = [t) (A20) 



Dchning the LLL bi-spinor as 

,/,, , , — I 

and using (|A16|) . Eq. (|A15[) can be rewritten in the compact form 

Co = ! ^ lll (p) [7 • P ~ E°] i> LLL (p) (A21) 

where = (po^Ps)- From the Lagrangian (|A21[) we obtain the Dirac equation for the LLL fermions 

[p-^-E°]^ LLL = 0, (A22) 

The Hamiltonian associated to £q is 

H = E Q +p 3 j 1 (A23) 

which is the free (1+1 )-D Thirring model Hamiltonian (2f| with the replacement of the fermion mass m by the rest- 
energy E°. Thus, in the chiral-condensate phase the LLL fermions behave as free particles in a reduced (1 + 1)-D 
space with rest energy proportional to the inverse magnetic length = \feH (see Eq. (|71[) ). 
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